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In the absence of decay, the conditional dynamics for an open system is often describable by a 
non-Hermitian Hamiltonian. This study investigates spin squeezing (SS) in non-Hermitian one-axis 
twisting (OAT) model. Somewhat surprisingly, SS close to the limit of Hermitian two-axis counter 
twisting (TACT) Hamiltonian is achievable for some parameters, which signihcantly improves upon 
the optimal value realizable by Hermitian OAT model. The drawback is like with all conditional 
schemes, it takes on average longer time to evolve into steady state, and the probability of no decay 
or success decreases as number of atoms (spins) increases. The result above for steady state SS in 
non-Hermitian OAT Hamiltonian is thus limited to small systems. For other parameter regimes, 
however, desirable SS arrives dynamically before steady state is achieved, with greatly shortened 
evolution time and enhanced probability of success, while still remain significantly improved over 
the limit of Hermitian OAT. 


I. INTRODUCTION 

Squeezed spin state (SSS) [1-3] is a symmetric state 
of spin-1/2 particles, whose fluctuation in one collective 
spin component perpendicular to the mean spin direction 
is smaller than the classical limit set by the summed fluc¬ 
tuations from independent spins pointing along the same 
direction. This reduced fluctuation manifests quantum 
correlations among individual spins. SSS has attracted 
considerable attention because of its potential applica¬ 
tion to improve the precision of quantum measurements 
[2-8] and to detect quantum entanglement [9-12]. Ac¬ 
cording to Kitagawa and Ueda [1], SSS can be dynam¬ 
ically generated from the product initial state with all 
spins up (or down) under the one-axis twisting (OAT) 
Hamiltonian Hqat = xJx ^^e two-axis counter 
twisting (TACT) Hamiltonian Htact = x{Jx ~ Jy)- 
The collective spin J (= with h = 1 here¬ 

inafter) is dehned in terms of the Pauli operator of the 
fc-th spin or pseudospin Ok, and y denotes the coupling 
strength between two spins. When SSS is applied to 
quantum metrology [2, 3], the property of interest is the 
squeezing parameter = iV(AJi)^/j(J)p < 1, where 
(AJj_)^ = {J\) — denotes the minimal fluctua¬ 

tion of a spin component perpendicular to the mean spin 
(J) = ((J.),(J,),(J.)). 

The optimal spin squeezing (SS) realizable theoret¬ 
ically for the OAT model is [1], which is 

reached at the time A better squeezing of 4/A 

[1] approaching the Heisenberg limit 1/A is achievable 
through TACT at a shorter time ln(4A)/(2A), which 
generally helps to mitigate accumulative influences from 
detrimental effects induced by particle losses and phase 
dephasing. Besides, unlike the situation encountered in 
OAT, the direction of optimal squeezing from TACT re¬ 
mains fixed during time evolution and is independent of 
system size (number of spins A) [13]. Despite of its bet¬ 


ter performance in SS, two body interactions capable of 
facilitating TACT do not occur naturally in most systems 
of interest. Although many proposals have been put for¬ 
ward to implement TACT models [13-18], no experimen¬ 
tal realizations have been reported so far. In contrast, SS 
from OAT has been proposed and demonstrated in vari¬ 
ous systems [4, 5, 7, 19-26]. Therefore, many subsequent 
studies aimed at improved SS in general OAT models are 
proposed [13, 27-30]. 

A recent study [31] reports the surprising finding 
that the non-Hermitian TACT model can realize slightly 
stronger SS than its Hermitian counterpart, which is 
counter intuitive as damping is always viewed as causing 
damage to quantum coherent processes like SS. Inspired 
by the desire for a better understanding of SS in non- 
Hermitian models, we carried out this investigation of 
the non-Hermitian OAT model. We find surprisingly the 
non-Hermitian OAT model may be even more favorable. 
In addition to being more readily realizable experimen¬ 
tally, it provides SS approaching the limit of TACT model 
l/^i Heisenberg limit) or even slightly better with 
a smaller coefficient when the system parameter is op¬ 
timal, which is significantly better than the optimal SS 
available from Hermitian OAT. As will be shown later in 
this study, near optimal squeezing, e.g., with scaling like 
^2 ^ can be reached within significantly shorter 

time at other system parameter values. Combining our 
results with that of Ref. [31], we come to the conclusion 
that the presence of dissipation can indeed improve the 
degree of squeezing, independent of the mechanism used 
to produce squeezing (OAT or TACT). 

This paper presents our investigation of enhanced SS 
in the non-Hermitian OAT model. It is organized as 
follows. Following this introduction section, the next 
two sections respectively discuss steady state SS and dy¬ 
namically generated optimal squeezed states in the non- 
Hermitian OAT model. We compare the above results to 
the corresponding ones from the OAT and TACT models. 
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We conclude in the last section. 


II. STEADY STATE 

Our model is built on the collective OAT interaction 
of an ensemble of N two-state atoms (i.e., a collection 
of pseudo-spin 1/2 particles) with up and down states 
{| t)j| i)}- Assuming a finite lifetime (I/ 7 ) for atoms 
in state | t); the density matrix p for atoms satisfies the 
master equation: dp/dt = —i[HoxT, p] + where the 
super-operator 

= ( 1 ) 

k 

with the decay rate 7 and the Pauli operators of the 
/c-th atom As in the non-Hermitian TACT model 
of Lee et al. [31], a finite and tunable value for 7 can 
be engineered through coupling of state | t) to an un¬ 
stable auxiliary state |a). The Pauli operators in the 
above Eq. (1) are ct/" = (Cfc )'*' = | t)kk{a\, hence 
^k^k ~ \ t)fcfc(t I- Conditioned on the absence of a de¬ 
cay event [31], whose probability is given by P = 

[32] for independent atomic decay as modeled here, one 
can remove the “real” decay term p<J^ in Eq. (1) and 
obtain jC{p) = -"f{N-fp +pN-f)/2, where Ni = 10fcfc{/| 

for the state label I = {t:i} denotes the atom-number 
operator. The conditional master equation now becomes 
dp/dt = —i{HesP — with the effective Hamilto¬ 

nian 

i/eff = xJl - *7iVt/2, (2) 

with N/2. 

The Hermitian OAT term can be realized for instance 
in trapped ions [20, 33] or cavity QED [34], starting 
with a coupled Hamiltonian of the form % = Aa^a -I- 
g {a -\- a^) Jx, where a and denote annihilation and 
creation operators of phonons in trapped ions or photons 
in cavity QED. For a trapped ion system, this is real¬ 
ized through the two photon interaction of ions with two 
lasers of opposite detunings [20, 33]. In cavity QED, it 
is implemented by the interaction of a single cavity pho¬ 
ton mode with atoms driven by a pair of coherent laser 
fields [34]. The evolution governed by H is described 
by propagator U {t) = g-i-f (*)<>■ Jx g-if* ^ 

with f (t) = i (g/A) — 1), and X(t) = 

- (g^/A) [t + i - 1) / A] [35]. In the weak cou¬ 
pling limit A 3 > g, / (t) is negligible, so we arrive at 
an effective Hamiltonian Hgs = —g^Jx/^- Beyond the 
weak coupling regime, ions (atoms) are strongly entan¬ 
gled with the vibrational motion (cavity photons). At 
times t = 2 /c 7 r/A, the vibrational motion (cavity mode) 
returns to its original state, and the propagator reduces 
to U (t) = 

The presence of the non-Hermitian term causes all 
eigenvalues of the model Hamiltonian (2) to possess neg¬ 
ative imaginary parts. The state with the largest imag¬ 
inary part becomes the steady state of the system as 








(d) 


FIG. 1. (Color online)(a) Spin squeezing parameter for 
the steady state of our non-Hermitian OAT model Hamilto¬ 
nian (2) as a function of 'y/x- The dashed (red) and dash- 
dotted (black) lines denote the optimal squeezing parameter 
from OAT and TACT respectively; The blue disk marks the 
parameter 7 /x = 0.4673, where the optimal squeezing is re¬ 
alized. The black diamond marks the parameter 7 /x = 1.19, 
where the squeezing limit of TACT is achieved. The red 
square marks the parameter 'y/x = 1.6393, where the de¬ 
gree of squeezing lies between the theoretical limits of the 
OAT and TACT, (b) Time evolution of the squeezing param¬ 
eter at the corresponding 7/7 values labeled by markers 
in Fig. 1(a), starting from the initial spin state jf-f • • • 4)- 
The dashed (red) line denotes the corresponding result for 
the OAT model; The blue, black, red solid lines denote the 
results for parameter 7 /x = 0.4673,1.19, 1.6393, respectively, 
(c) The corresponding probability of success P{= as 

a function of evolution time, (d) Quasi-probability distribu¬ 
tion Q {0, ip) in the Jx-Jy plane for the three states labeled as 
(i)-(iii) in (b). 7 /x = 1.19. Jz is chosen to be along the direc¬ 
tion pointing into the page with 8 and (p the corresponding 
polar and azimuthal angles. Here the illustrative calculations 
are carried out for N = 20. 


it eventually becomes the lone survivor after a suffi¬ 
cient time of evolution. We now compute its squeezing 
properties . The Hamiltonian (2) maintains parity sym¬ 
metry in Jx and Jy, which assures {Jx) = (Jy) = 0. 
The mean spin therefore points along the z-axis. The 
squeezing parameter for its steady state is determined by 
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with J_L = Ja^cosa + JySina 
lie in the x-y plane. Given the system size N, the ratio 
7 /x is the only parameter governing the properties of the 
Hamiltonian (2). Depending on its value, the squeezing 
parameter exhibits different behavior due to the com¬ 
petition between the OAT interaction and the dissipa¬ 
tion. As mentioned earlier, 7 is tunable [31] and can 
assume whatever A-dependence needed. For essentially 
all atomic systems, spontaneous emission can be added 
to an atomic state from induced off-resonant coupling to 
an unstable state. 

Figure 1(a) presents the squeezing parameter [blue 
solid line] obtained numerically as a function of "f/x- 
The actual calculation is carried out using exact diag- 
onalization of the non-Hermitian Hamiltonian (2). Over 
a wide range of 7 /% [shaded area], the steady state of 
the non-Hermitian Hamiltonian (2) is found to possess 
more SS than the optimal SS afforded by the Hermi- 
tian OAT model [red dashed line]. In some parameter 
regime, the degree of SS is found to even surpass the 
much improved limit provided by the Hermitian TACT 
model [black dash-dotted line]. 

To further confirm their validity, we simulate the dy¬ 
namical non-Hermitian evolution governed by Hamilto¬ 
nian ( 2 ) for the squeezing parameter with 7 /x set at 
the marked values in Fig. 1(a) for the same initial state 
with all spins down [),). The corresponding results are 
shown in Fig. 1(b). The detailed dynamics depend on 
the initial state, but the steady state squeezing proper¬ 
ties do not. The initial state with all spin down is the 
only stable N atom state, which represents the natural 
starting point for the conditional dynamics described by 
the non-Hermitian Hamiltonian (2). At short times, the 
non-Hermitian term has little effect and the squeezing 
parameter is observed to decrease according to OAT. 
After reaching the optimal point of minimum squeezing, 
which is essentially equal to limit of Hermitian OAT, its 
value starts to increase as the non-Hermitian term comes 
into play. This clearly shows up as the two curves for Her¬ 
mitian [red dashed line] and non-Hermitian OAT [solid 
lines] start to deviate from each other. The squeezing 
parameter for non-Hermitian OAT continues to de¬ 
crease as time goes on and finally approaches its steady 
state value. The degree of optimal squeezing and the 
time to achieve it depends on the parameter 7 /%. For 
7 /x = 0.4673 [blue disk in Fig. 1(a)], the squeezing pa¬ 
rameter [blue solid line in Fig. 1(b)] reaches its opti¬ 
mal value at tmin 2 /x, and remains steady at this level 
as time goes on, although the likelihood for a decay event 
to destroy coherence increases. For a larger 7 /x = 1.19 
[black diamond in Fig. 1(a)], a weaker squeezing equiv¬ 
alent to the TACT limit [black solid line in Fig. 1(b)] 
is obtained at a shorter time tmin 1/x- Different from 
the above two cases, where the time for the squeezing 
parameter to reach its steady value G matches with tmin, 
for 7 /x = 1.6393 [red square in Fig. 1(a)], the optimal 
value of [red solid line in Fig. 1(b)] arrives at an ear¬ 
lier time tniin 0.3/x before the system settles down to 





FIG. 2. (Color online) The A-dependence of the various SS 
quantities for the steady state, (a) squeezing parameter 
(b) the angle a?min between the optimal squeezing direction 
with respect to a;-axis; (c) time to reach steady state x^s; (d) 
total evolution time xts/F; all for 7 /x = 1/(0.03A). The 
red dashed and black dash-dotted lines denote respectively 
the corresponding results of the Hermitian OAT and TACT 
models. 


the steady state, i.e., G > ^min- 

The SS we study in this work for the non-Hermitian 
OAT model is only experimentally accessible when no de¬ 
cay event occurs, thus the probability of success P also 
represents an important consideration. Figure 1(c) shows 
the probability of success P as a function of evolution 
time. By comparing the results from the three different 
7 /x values, we find that the degree of squeezing decreases 
as 7 /x increases, while the corresponding P increases. 
For instance, at 7 /x = 1.19, the probability of success 
P when the optimal squeezing is achieved is about 0.4, 
i.e., on average for two out of five experimental runs, no 
decay occurs before reaching optimal squeezing; While at 
7 /x = 1.6393, the corresponding P ~ 0.6. These rates 
of success imply that the non-Hermitian scheme is fea¬ 
sible, certainly for the small size of A = 20 numerically 
evaluated here. 

The squeezing process can be intuitively illus¬ 
trated by the evolution of quasi-probability dis¬ 
tribution Q{9,ip) for the state \tp{t)), which is 
determined by its projection onto the coherent 
spin state J^, (/?)=(cos(0/2)j f) -I- e*‘^sin(0/2)j i))^^ , i.e., 
Q{9,ip) = \{6,ip\tjj{t))\'^, as shown in Fig. 1(d). Shearing 
of the initial isotropic uncertainty distribution results in 
reduced spin variance along one direction. Defining amin 
as the angle between the optimal squeezing direction with 
respect to x-axis, it is seen that Omin tends zero under 
the influence of the non-Hermitian term. 

To check if the non-Hermitian scheme for improved SS 
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remains valid at larger atom numbers we carried out 
calculations which clearly reveal the scaling of various SS 
quantities with respect to N as shown in Fig. 2. The pa¬ 
rameter 7 /x is set at 1/(0.03iV), and the results are fitted 
to give a degree of squeezing approaching the TACT limit 
4/N, as depicted in Fig. 2(a). The squeezing direction 
is almost fixed along x-axis or amin — 0 as shown in Fig. 
2(b), representing a significant advantage over the Her- 
mitian OAT model. The drawback is that the optimal 
squeezing time tmin) which is essentially the same time 
for reaching steady state tg, is around 0 . 1 /% ~ 1 /x and 
it shows a weak dependence on N, slightly increases in 
the range of A^ = 100 to A^ = 10000, as illustrated in Fig. 
2(c). This is in contrast to the Hermitian OAT model, 
whose squeezing time decreases as N increases. In practi¬ 
cal implementations, this could present a serious obstacle 
for systems with short coherence times. The probability 
of success P decays exponentially with TV, like other con¬ 
ditional schemes. For small TV < 1000, P is found to 
decay slightly more rapidly than for TV > 1000. As our 
model is only experimentally accessible when no decay 
event occurs, a low success rate implies more measure¬ 
ment runs, which directly translates into longer times. 
To reach a success rate P, the total evolution time for 
success will have to be around xts/P- For instance, for 
P = 0.2, five experiments on average must be conducted 
to give a successful event, thus, the total evolution time 
is Syts = xts/P, with each experiment cost time xts- 
As shown in Fig. 2(d), the total evolution time for suc¬ 
cess xts/P increases as atom number TV increases. For 
TV > 1000, it becomes larger than 10. 


III. OPTIMAL SQUEEZED STATE 

In previous discussions, we focus on the steady state at 
a specific parameter of 7 /x = 1/(0.03TV), which provides 
a significantly improved SS with the degree of squeezing 
equal to the Hermitian TACT model. The price to pay is 
the prolonged evolution time to settle into steady state 
and the low probability of success which further decreases 
as TV increases. To overcome these problems, we search 
for other parameter regimes. 

For instance, we show squeezing parameter obtained 
from non-Hermitian evolution at parameters 7/x = 0.1 
[blue solid line] and 7/x = 0.5 [blue dotted line] for 
TV = 1000 in Fig. 3(a). For both cases, we find their 
reach the minimal values [marked by arrows] before 
settling down to the steady values, i.e., tmin < ts- Al¬ 
though the optimal squeezing at these two values for 
the parameter 7/x does not surpass the squeezing at 
7 /x = 1/(0.03TV), which approaches the TACT limit, the 
benefits for these two cases reside in their much short¬ 
ened evolution times. To learn the dependence of the 
optimal squeezing parameter on the parameter 7/x, 
we show the relations between them in Fig. 3(b). In 
a wide parameter range, the degree of squeezing for our 
non-Hermitian OAT model again is found to surpass that 




7/x 



7/x 


FIG. 3. (Color online) (a) The squeezing parameter ^ ob¬ 
tained from non-Hermitian evolution at 7 /x = 0.1 (blue solid 
line) and 7 /x = 0.5 (blue dotted line), starting from initial 
state IF! • • • 4 )- (b)-(d) denote respectively the scaling of var¬ 
ious SS quantities with y/x- (b) optimal squeezing parameter 
Cmin; (c) optimal squeezing time x^min; and (d) total evolu¬ 
tion time xtmin/P- Filled squares and disks denote the results 
from non-Hermitian OAT model at 7 /x = 0.1 and 7 /x = 0.5, 
respectively. The red dashed and black dash-dotted lines de¬ 
note the corresponding results from OAT and TACT. All cal¬ 
culations are for N = 1000. 


of the Hermitian OAT model, which implies that the per¬ 
formance of our non-Hermitian scheme is insensitive to 
noise induced parameter fluctuations. With smaller 7/x, 
we find the degree of squeezing improves, again at the 
cost of longer evolution times, whether for a single run 
[Fig. 3(c)] or for the total time needed for success [Fig. 
3(d)]. 

The results above therefore impose a trade-off between 
the degree of squeezing with evolution time. Reaching 
the optimal degree of squeezing requires a delicate bal¬ 
ance between the two. To demonstrate this more clearly, 
we compare the scaling of various quantities with TV 
at two values of parameter 7/x = 0.1 and 0.5 as pre¬ 
sented in Fig. 4. For 7/x = 0.1 [blue squares], we find 
the squeezing parameter scales as 1.2TV“^/® [Fig. 4(a), 
blue solid line], the corresponding evolution time for 
a single run xTmin scales approximately as 17.4TV“^/^ 
[Fig. 4(b), blue solid line], and the total evolution time 
to success xtmin/P decreases as TV increases [Fig. 4(c)]. 
For the larger 7/x = 0.5 [blue disks], the degree of 
squeezing becomes less with scaling 1.2TV“^/^ [Fig. 
4(a), blue dashed line], the evolution time both for a 
single run 5.4TV“^/^ [Fig. 4(b), blue dashed line] and 
for the total [Fig. 4(c)] become shorter. Specifically, for 
TV = 10^ spins, we obtain a degree of squeezing equal 
to lOlog^o^^ = —31.1dB for the former at an evolu¬ 
tion time of xTmin = 0.0354, while for the latter case, 
—29.0dB squeezing is reached at x^min = 0.0112. They 
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FIG. 4. (Color online) The dependence of various SS quan¬ 
tities on N. (a) Optimal squeezing parameter ^mini (b) op¬ 
timal squeezing time (c) total evolution time required 

for success xtmin/P- Dotted markers denote the results from 
non-Hermitian OAT model with 7 /x = 0.1 (blue squares) and 
0.5 (blue disks). The red dashed and black dash-dotted lines 
denote the corresponding results from Hermitian OAT and 
TACT models, respectively. 


can be compared to the more standard results of —26.2dB 
at Xtmin = 0.00254 for the Hermitian OAT model, and 
—34.1dB at xtmin = 0.000445 for the Hermitian TACT 
model. 


IV. CONCLUSION 

In conclusion, we investigate SS in non-Hermitian OAT 
model, which is conditional on the absence of decay event. 
The ratio of the dissipation rate 7 to the OAT interac¬ 
tion strength y, 7 /x, is the only tunable parameter in 
this model, whose value determines the complete squeez¬ 
ing behavior. At 7 /x = 1/(0.03A^), we find steady state 
SS can reach the squeezing limit of TACT, which repre¬ 
sents a significant improvement over the Hermitian OAT 
model. This is achieved at the expense of prolonged evo¬ 
lution tiirre and low probability of success P conditional 
on no decay event. The former iircreases slightly with V, 
while the later decays exponentially with V, like other 
conditional schemes. This combines to give an overall un¬ 
favorable scaling with TV. Thus while encouraging, our 
scheme is perhaps only applicable to systems with small 
N. Furthermore, by investigating the non-Hermitian dy¬ 
namics at other parameter regimes where optimal SSSs 
arrive before steady states, we find optimal squeezing 
time can be greatly shortened and success rate enhanced, 
at the same time impressive degrees of SS significantly 
beating the OAT model are maintained. Our work high¬ 
lights potentially fruitful applications of non-Hermitian 
OAT to small samples of coupled spins, including multi¬ 
ple component atomic condensates where non-Hermitian 
OAT model Hamiltonian can be engineered and experi¬ 
mentally implemented. 
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